Star-quantization of an infinite wall 
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Abstract: In deformation quantization (a.k.a. the Wigner-Weyl-Moyal formulation of 
quantum mechanics), we consider a single quantum particle moving freely in one dimension, 
except for the presence of one infinite potential wall. Dias and Prata pointed out that, 
surprisingly, its stationary-state Wigner function does not obey the naive equation of motion, 
i.e. the naive stargenvalue (*-genvalue) equation. We review our recent work on this problem, 
that treats the infinite wall as the limit of a Liouville potential. Also included are some new 
results: (i) we show explicitly that the Wigner- Weyl transform of the usual density matrix 
is the physical solution, (ii) we prove that an effective-mass treatment of the problem is 
equivalent to the Liouville one, and (iii) we point out that self-adjointness of the operator 
Hamiltonian requires a boundary potential, but one different from that proposed by Dias and 
Prata. 

PACS Nos.: 03.65.-w, 03.65.Ca, 03.65.Ge. 

Resume : French version of abstract (supplied by CJP if necessary) 
[Traduit par la redaction] 

1. Introduction 

Deformation quantization is also known as the Wigner-Weyl-Moyal formulation of quantum me- 
chanics, and by other names. 1 We'll refer to it here as star-quantization. In it, as in classical mechanics, 
observables are realized as ordinary functions (and distributions) on phase space. They are multiplied, 
however, using a non-commutative, associative product, known as the star product (^-product). 

Star-quantization can be understood as the Wigner- Weyl transform of the quantum dynamics of the 
density matrix. On the other hand, it is an autonomous way of doing quantum mechanics. Quantum 
systems should therefore be treatable in it, without reference to operators and wave functions, or to path 
integrals, or to any other quantum formulation. That is, physical systems should be star-quantizable. 

This note is concerned with the star-quantization of a very simple system, one particle moving 
freely along the negative x-axis, in the presence of an infinite potential wall at x = 0. Naive star- 
quantization does not work for it [2]. We would like to understand why it fails, and how the problem 
can be fixed. 

Our considerations will be restricted to stationary states of the system. As emphasized in [3], where 
several examples were worked out, such stationary states are described by time-independent Wigner 
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functions that can be found as solutions of the so-called *-genvalue equations. 

We first introduce the *-genvalue equations in section 2, and solve them for the Wigner function of 
a free particle moving on the entire ir-axis. In the following section, the original problem [2], occurring 
when an infinite barrier is present at x = 0, is described. In section 4, we review our treatment [4] 
of the star-quantization of an infinite wall as the limit of an exponential potential. An alternative to 
the *-genvalue equation was found that can be treated in the naive way: it can be solved first, and 
then the boundary conditions can be imposed later, to yield the physical Wigner function. In section 5, 
we note that the system can be treated using an effective mass. Interestingly, however, we show that 
the corresponding *-genvalue equations are equivalent to those for the exponential potential. The ad 
hoc resolution of the problem proposed by Dias and Prata [2] (see also [5, 6]) is then described, in 
section 6. We have found in operator quantum mechanics a motivation for a result similar in form to, 
but in conflict with theirs. This is sketched in the same section 6. In [7], we plan to work out further 
consequences in star-quantization of this observation. 

2. Free particle: *-genvalue equation & Wigner-Weyl transform 

Consider a single particle moving on the x-axis, so that its phase space is R 2 , with coordinates 
x,p. In star-quantization, the Wigner function encodes all observable information about the quantum 
state. Stationary states are described by time-independent Wigner functions p(x,p), that can be found 
as solutions to the so-called *-genvalue equations 

H*p = p*H = Ep . (1) 

Here, H = H(x,p) is the classical Hamiltonian, E is the energy, and the Gronewold-Moyal ^-product 
is defined by 

f{x,p)*g(x,p) = ex.p[ih(d x d p >-dpd x <)/2}f(x,p)g(x',p')\ {x , p , )={xp) . (2) 

Real solutions of these equations will describe quantum states, but they will be mixed states, in gen- 
eral. In this work, we will only consider pure states, whose Wigner functions obey the additional 
★-projection condition 

Pi*Pj = 6 lj p j /h, or p a *pj3 oc S(a- (3) pp/h , (3) 

for states labelled by discrete parameters i,j,..., or for non-normalizable states labelled by continuous 
parameters a,f3,.. ., respectively. 

To illustrate, let us first consider a free particle (see the Appendix of [4]). The Hamiltonian is 
H = p 2 , if we set 2m = 1 for simplicity. The *-genvalue equations (1) then reduce to 

(p-ihd x /2) 2 p(x,p) = k 2 p(x,p) , (4) 

where we have put E = k 2 , and we will work henceforth with Ti = 1. The imaginary part of this 
equation ispd x p = 0, so that when p ^ 0, d x p = 0, but not when p = 0. This leads to the ansatz 

p(x,p) = f(p) + S(p)g(x). (5) 

In the real part of (4), (p 2 - 9^/4) p = k 2 p, it yields 

(P 2 - k 2 ) f(p) - S(p)(k 2 + d 2 /4) g(x) = . (6) 
Considering p — gives 

g(x) = benp(2ikx) + b* exp(— 2ikx) , (7) 
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where p = p* has been imposed. Then (6) reduces to (p 2 — k 2 ) f(p) — 0, solved by 

/(p) = a+6(p-k) + a.5{p + k), (8) 

with a± arbitrary real constants. The general result is therefore 

p = S(p) { bexp(2ikx) + b* exp(~2ikx) } + a + 5(p — k) + aS(p + k) . (9) 

To restrict to pure-state Wigner functions, we impose p* p oc 5(0) p, a special case of (3). We find 
the constraint 

\b\ 2 = a+a_ => b = y/a^aZe^ , <p e M . (10) 
The general pure-state solution to the free-particle *-genvalue equation (4) is therefore 

p = a + S(p — k) + ci-5(p + k) + 2^/a + a_ 5(p) cos (2kx + <f) . (11) 

On the other hand, the Wigner function can be calculated as the (normalized) Wigner- Weyl trans- 
form of the density matrix built from the known wave functions: 

/#] := - / dye- 2 ^yj(x + y)r^~y) ■ (12) 

We denote it as p[tp] here to emphasize that it is calculated from known input wave functions, and so 
is not a result of working within the star-quantization formulation. With the pure-state wave function 
■0 = a + e tkx + a.-e~ lkx , eqn. (12) yields 

p[V>] = \a+\ 2 5(p-k) + \a-\ 2 5(p + k) + 5{p){a* + a^e- 2lkx + a+a*_e 2lkx } . (13) 
Comparing (13) with (11) reveals a one-to-one correspondence. All is well for the free particle. 

3. Infinite wall: *-genvalue equation vs Wigner-Weyl transform 

Now consider the Hamiltonian H = p 2 + V(x), with potential energy 

V{x) = { °' X< °> (14) 
w [ oo , x > . 

We first work out the Wigner-Weyl transform of the known Schrodinger wave functions 

tp(x) = 9(-x) sin(kx) . (15) 

Here 6(x) denotes the Heaviside step function. From (12) we find 

p[ip] oc 0(-x)p(x,p) , (16) 
with 

This can be obtained from a special case of the free-particle Wigner function (1 1) by the replacement 
S(p — Po) — * sin[2x(p — p )]/[2(p — p )]. The latter vanishes at x = 0, so that p{0,p) — 0, and 
becomes the former when x — > — oo, so that a free-particle result is found far from the infinite wall. 

We will show in the following section that (16, 17) describe the physical Wigner function, the one 
that should be found in the star-quantization of an infinite wall. 

To solve the *-genvalue equations (1) with the Hamiltonian determined by (14), we can try to follow 
the Schrodinger treatment of this system, by restricting to x < 0, and then imposing the boundary 
condition p(0,p) = 0. For x < 0, the *-genvalue equation is that of a free particle, however, eqn. (4). 
The problem is that the Wigner function given by (16, 17) does not satisfy this equation. 
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4. Infinite barrier as the limit of an exponential potential 

To study this problem, we treated the Hamiltonian of the previous section as the a — > oo limit of 

H a =p 2 + e 2ax . (18) 

Star-quantization of this system had already been carried out at a = 1, in [3]. The a-dependence was 
easily reinstated. The physical Wigner function is given by (cf. eqn. (98) of [3]) 

Pa oc f c ds [^gr] r(^(p-fc)-s)r(^(p + fc)-s) 

xr(g(p-fc)- s )r(g( P + fc)- s ) . (19) 

The contour C in the s plane runs from — ico to +ioo, just to the left of the four poles on the imaginary 
s axis at ± ^{p ±' k). The right-hand side of (19) is the Mellin-Barnes type integral definition of the 
Meijer G-function (see §5.3. of [8]), so 



p a oc G 



{2af 



i(p — k) i(p + k) —i{p — k) —i(p + k) 



2a ' 2a ' 2a 2a 



(20) 



First, consider the a — > oo limit of the Wigner function. For x > 0, w — e 4ax /(2a) 4 — > oo, but 
GaUM ) ~~ * exponentially, according to the asymptotics described in §5.4.1 of [8]. 

For x < 0, the contour C can be closed such that all 4 poles on the imaginary s-axis, at s — 
i(±p ±' k)/2a, are surrounded. First note that 

lim ( 2a )- 24 ( ± P ± ' fe )/ Q = 1 . (21) 

We can now apply the residue theorem and use T(z) = T(l + z)/z, to find p a (x,p) goes as 

e iix( P -k) /(-2ik)2(ip- ik){2ip) + e 2lx (P+ k y(2ik){2ip)2(ip + ik) 
+ e -2ix(p-k)/2(-ip + ik)(-2ip)(2ik) + e- 2tx (P+ k y(-2ip)2(-ip- ik)(-2ik) , (22) 

as a — > oo, up to a factor depending on a, but independent of x. But this is proportional to 

2ix(p-k) ( 1 ]_\ i 2ix{p+k) ( i i A 

e \kp k(p-k)) ^ e \kp k(p+k)J 

_ -2ix(p-k) ( 1 1_\ _ -2ix(p+k) ( 1 ]_\ (2 T.\ 

e \kp k(p-k) J \kp k(p+k) J ■ v"-* 

Comparing with (17), we see that for x < 0, we have shown that linia^oo p a oc p. 

Combining the results for x < and x > 0, we have that linio^oo p a oc 6(— x) p. That is, the limit 
of the Wigner function for the exponential potential is 9(— x) p, showing that (16) and (17) describe the 
physical Wigner function for the infinite potential wall. The authors of [2], e.g., made this assumption. 

Our goal is to find a dynamical equation that this physical Wigner function satisfies, instead of the 
★-genvalue equation, which it does not. The new equation [4] can then be used to solve for the physical 
Wigner function, i.e., the new equation can be used to carry out a true star-quantization of the infinite 
potential wall. 

Following [3], we see that (18) yields a *-genvalue equation H a * p a = k 2 p a with imaginary and 
real parts 

e -2ax d x p a (x,p) = - ^[p a (x,p + ia) - p a (x,p- ia)} , 
e- 2ax (p 2 -k 2 -\dl)p a {x,p) + \ [ Pa {x,p + ia)+p a {x,p-ia)] = 0, (24) 
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respectively. These two equations can be recombined into a difference equation: 2 

0= [f-^) Poi {x,p) + l(^-) 



2 

p a (x,p+2ia) — p a (x.p) | p a (x.p— 2ia) — p a (x,p) 
p+ia p—ia 



ip 



[p a (x,p + ia) - p a {x,p - ia)] + [p a (x,p + ia) + p a (x,p - ia)] . (25) 



In [4], we showed that Pa (x,p ± ia), p a (x,p±2ia) could be traded for <9™ p a (x,p), n= 1,2,3,4, 
to arrive at a differential equation, with a simple a — > oo limit, 

d 4 xP (x,p)/16 + (p 2 + k 2 )d 2 x p(x,p)/2 + (p 4 -2k 2 p + k 4 )p(x,p) = 0, (26) 
valid for x < 0. Recalling that E =: k 2 , the new equation can be rewritten as 

(p 2 - E)* p*(p 2 - E) = . (27) 

It is easily verified that p of (17) satisfies the new equation. In [4] it was also shown that the 
new equation applies to other systems built from infinite walls or wells: the infinite square well and 
the delta-function well. Furthermore, it was shown how an additional, non-singular potential could be 
included in the framework. 

The proposal made in [4] was to star-quantize the infinite-wall system by solving the new equation, 
and then imposing the boundary conditions, much as for solutions of the Schrodinger equation. Some 
progress in this direction was made in [6]. 

5. Effective mass treatment 

Let us mention another possible approach. It is clear that the free Hamiltonian H = p 2 encodes 
nothing of the dynamics of the infinite potential wall at x = 0. In the Schrodinger treatment, imposing 
the boundary conditions suffices to yield the required physics, but not in star-quantization. Can the 
dynamics of the wall be described by a modified Hamiltonian? 

One could try to incorporate the dynamics in the kinematics, by assigning to the particle an effective 
mass that blows up for x > 0. In the spirit of the last section, such a Hamiltonian can also be considered 
as a limit: 

H' eS := 9(-x)p 2 = lim H' eS a := lim (l + e 2ttx ) _1 p 2 . (28) 
We will now show that 

H' cSi a (x' , p) * p' cS (x',p) = k 2 p' cS (x',p) (29) 

is equivalent to the *-genvalue equation with Hamiltonian H a of (18), already considered. 
Using (2), one can show that (29) reduces to 

(p- t d x/ /2) 2 p>^ a (x',p) = k 2 (l + e 2ax ' + ^)p'^ a (x' 7 p). (30) 
Writings' = x + log(-fc 2 )/(2a) yields 

(p-id x /2) 2 p a (x,p) - (k 2 -e 2ax+lad *)p a {x,p) , (31) 

however, where we put p a (x,p) :— p' cS a (x',p). The transformed equationis just H a (x,p)*p a (x,p) = 
k 2 p a (x,p), the *-genvalue equation for (18). Since x' — x — > in the a — > oo limit, (29) must lead to 
the same eqn. (27) we already found above. 

2 At a = 1, this is a slight correction of eqn. (104) of [3]. 
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6. Self-adjoint extensions and boundary potentials 

Let us reconsider the operator treatment of the system, to try to understand why pure star-quantization 
fails. Perhaps the problem will turn out to be the Wigner-Weyl transform of one already present in the 
operator formulation. 

A self-adjoint Hamiltonian H will have a spectral decomposition, H = "Y^ E , \E')(E'\, roughly 
speaking. Then for a pure-state density matrix p = \E) (E\ built from the stationary energy eigenstate 
\E), we necessarily have Hp = pH — Ep. If the Wigner-Weyl transform W works, these become 
the *-genvalue equations (1). 

The key realization is that the free Hamiltonian operator H := p 2 is not self-adjoint on the negative 
x-axis, even if it is Hermitian (see [9], e.g., and references therein). We believe that this is the root 
cause of at least part of the problem encountered with naive star-quantization of the infinite wall. 

The free Hamiltonian H does have self-adjoint extensions, however, as can be determined by cal- 
culating its von Neumann deficiency indices. The extensions amount to including a point interaction at 
x = 0. The role of the point interaction is to enforce a (Robin) boundary condition on the Schrodinger 
wave function, involving a real length L: 

V>(0) + LV'(O) = (V':=#Mr). (32) 

It is easy to see that such a potential takes the form Vl (x) = S(x) — LS'(x). 

A specific interaction of this type was prescribed in an ad hoc way by Dias and Prata, for the 
boundary condition ^>(0) = 0. It is interesting to note, however, that their prescription was for a 
boundary potential proportional to 5'_(x), a regularized version of 5'(x). Apparently, this contradicts 
the standard operator treatment, since the Dirichlet boundary condition corresponds to L = in (32), 
which relates to Vl=o{x) = 5(x). 

Clearly, the situation needs clarification. We hope to report further progress on this problem in [7]. 
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